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The  present  project  is  aimed  at  developing  validated  models  for  predictions  of  static  stability  of 
combustion  in  augment  or  systems,  and  further  at  the  design  of  novel  flame  holding  concepts.  In 
the  present  reporting  period,  we  performed  large-eddy  simulations  of  the  Stanford  vitiated  flow 
experiment.  We  also  proposed  and  computationally  tested  the  concept  of  Damkohler  number 
similarity  for  assessing  flame  stability  in  low  speed  experimental  facilities.  We  further  acquired  a 
V-gutter  geometry  and  the  geometry  of  the  new  AFRL  flame  holder  test  facility  and  performed 
cold  flow  and  reactive  simulations  of  a  specified  flow  condition  as  a  first  point  of  a  flame  stability 
map. 

The  experimental  work  has  focused  on  improving  the  performance  of  the  bluff-body  augmen¬ 
tor  system  in  preparation  for  detailed  combustion  studies.  Improvements  consist  of  modifications 
to  the  bluff-body  fuel  injection  system  and  the  implementation  of  a  turbulence  grid.  We  have 
also  investigated  the  injection  of  hydrogen  peroxide  (H2O2)  in  the  post  combustion  stream.  The 
use  of  a  grid  with  higher  turbulence  conditions  is  used  as  a  proxy  for  higher  speed  flow,  while  the 
use  of  hydrogen  peroxide  is  used  to  study  the  effect  of  OH  radical  concentration  on  flame  sta¬ 
bility.  The  work  described  here  has  been  conducted  by  post  doctoral  researcher  Wookyung  Kim 
with  assistance  from  graduate  students  Hyungrok  Do  and  Seong-kyun  Im  under  the  supervision 
of  Prof.  Godfrey  MungaL  Details  of  the  efforts  axe  provided  below. 
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The  overall  objective  of  the  project  is  to  develop  well  tested  and  validated  models  that  can  predict 
processeses  for  static  stability  in  angmentor  systems  leading  to  a  better  design  of  advanced, 
compact,  stable,  and  efficient  angmentors.  The  emphasis  of  the  proposed  phase  II  work  is  to 
study,  fundamentally  understand,  and  model  the  chemical  effect  of  vitiation  on  static  stability  in 
augmentor  systems  with  computational  and  experimental  methods,  and  to  enhance  and  validate 
the  predictive  capabilities  of  the  commercial  multi-physics,  unstructured  large-eddy  simulation 
code  CharLES  distributed  by  CASCADE  Technologies.  We  will  also  propose  and  investigate 
new  flame  holder  designs  with  improved  static  stability. 

The  project  focusses  mainly  on  three  different  fundamental  aspects: 

Task  1:  Investigation  of  the  effects  of  vitiated  oxidizer  flow  composition  on  flame  stabilization. 
Experiments  are  carried  out  in  the  Stanford  experimental  rig  to  study  flame  stabilization 
in  vitiated  air*  Different  additives  to  the  vitiated  flow  are  considered  to  study  the  ef¬ 
fects  of  several  species,  including  OH,  NO,  and  NO2,  in  more  detail.  The  experimentally 
investigated  flows  are  also  computed  with  LES* 

Task  2:  Computational  model  extension  to  consider  mixed- mode  augmentor  stabilization  and  model 
validation.  Model  extensions  that  allow  for  a  more  accurate  consideration  of  finite  rate 
chemistry  are  being  developed*  Specifically,  a  model  will  be  developed  that  can  predict 
auto-ignition  in  turbulent  environments,  partially  premixed  turbulent  combustion,  non- 
premixed  turbulent  combustion,  and  the  transition  between  the  three  regimes.  All  model 
developments  will  be  validated  using  the  data  from  our  own  experiments  and  from  the 
AFRL  HIT  experiments* 

Task  3:  Conceptual  design  of  new  augmentor  flame  holder  with  improved  static  stability.  We  will 
work  on  the  conceptual  design  of  new  augmentor  flame  holders  and  test  promising  concepts 
in  the  Stanford  experimental  facility  and  at  AFRL.  The  Damkohler  number  similarity  will 
be  demonstrated  first.  Pre-design  of  new  flame  holding  concepts  will  then  be  developed 
in  the  Stanford  vitiated  flow  experiment*  The  designs  will  finally  be  tested  and  further 
refined  in  the  AFRL  ASTAR  experimental  facility. 

There  is  essentially  no  change  from  the  originally  proposed  work*  However,  we  propose 
a  shift  of  emphasis  to  a  closer  collaboration  with  the  newly  planned  AFRL  experiments.  This 
requires  more  emphasis  on  simulations  of  the  AFRL  data  and  less  on  simulations  of  the  Stanford 
experiments  than  originally  planned.  We  have  now  a  more  concrete  plan  for  the  augmentor 
design  and  model  validation  that  was  worked  out  in  close  collaborations  with  Dr*  Kiel  from 
AFRL.  We  will  first  validate  the  models  with  the  Stanford  low  speed  rig  and  the  AFRL  high 
speed  rig,  specifically  using  the  new  AFRL  data  that  is  presently  being  acquired.  We  are 
presently  computing  stability  maps  of  the  AFRL  configuration,  for  which  experiments  will  be 
done  in  the  near  future.  The  newly  proposed  flame  holder  designs  will  also  be  assessed  in  close 
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collaboration  with  AFRL.  New  geometries  will  first  be  assessed  computationally  and  potential 
design  changes  will  be  performed.  The  design  will  then  be  tested  in  the  Stanford  low  speed 
rig  and  the  AFRL  high  speed  facility.  The  test  conditions  will  be  chosen  using  the  concept  of 
Damkohler  number  similarity  that  is  described  in  the  present  report.  This  will  test  both  the 
predictive  capability  of  the  models,  but  also  the  predictiveness  of  the  concept  of  Damkohler 
number  similarity  that  is  proposed  in  the  present  report. 

3  Status  of  Effort 

The  present  project  is  aimed  at  developing  validated  models  for  predictions  of  static  stability  of 
combustion  in  augmentor  systems,  and  further  at  the  design  of  novel  flame  holding  concepts.  In 
the  present  reporting  period,  we  performed  large-eddy  simulations  of  the  Stanford  vitiated  flow 
experiment.  We  also  proposed  and  computationally  tested  the  concept  of  Damkohler  number 
similarity  for  assessing  flame  stability  in  low  speed  experimental  facilities.  We  further  acquired  a 
V-gutter  geometry  and  the  geometry  of  the  new  AFRL  flame  holder  test  facility  and  performed 
cold  flow  and  reactive  simulations  of  a  specified  flow  condition  as  a  first  point  of  a  flame  stability 
map.  The  experimental  work  has  focused  on  improving  the  performance  of  the  bluff-body 
augmentor  system  in  preparation  for  detailed  combustion  studies.  Improvements  consist  of 
modifications  to  the  bluff-body  fuel  injection  system  and  the  implementation  of  a  turbulence 
grid.  We  have  also  investigated  the  injection  of  hydrogen  peroxide  (H2O2)  in  the  post  combustion 
stream.  The  use  of  a  grid  with  higher  turbulence  conditions  is  used  as  a  proxy  for  higher  speed 
flow,  while  the  use  of  hydrogen  peroxide  is  used  to  study  the  effect  of  OH  radical  concentration 
on  flame  stability. 

4  Accomplishments 

In  the  present  reporting  period,  we  have  worked  on  all  tasks  described  above.  In  the  experimental 
part,  we  have  made  significant  modifications  to  the  experimental  facility  to  make  the  flow  more 
homogeneous  and  more  relevant  to  realistic  augmentor  systems.  Data  has  been  acquired  in  the 
new  facility  for  validation  and  for  assessing  the  influence  of  individual  species  in  the  vitiated 
stream.  In  the  computational  and  modeling  part  of  the  project,  we  have  performed  validation 
simulations  of  the  Stanford  experiment.  We  have  then  used  this  case  to  propose  and  demonstrate 
the  concept  of  Damkohler  number  similarity.  Finally,  we  have  acquired  the  geometry  of  the  new 
AFRL  high  speed  test  rig,  and  we  have  performed  several  non-reactive  and  reactive  simulations 
to  compare  reactive  versus  non-reactive  flow  and  LES  versus  2D  RANS  simulations.  These 
accomplishments  are  presented  in  the  following. 
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4.1  Experimental  Study 

4.1.1  Modifications  to  the  Bluff-Body  Fuel  Injection  System 

The  overall  geometry  of  the  bluff  body  is  unchanged  from  that  of  Phase  I,  but  the  number  of 
fuel  injection  holes  on  its  sides  has  been  adjusted  from  two  centrally  located  holes  in  Phase 
I  to  ten  holes,  then  finally  settled  onto  a  six  hole  configuration.  The  primary  aims  of  these 
modifications  are  to  increase  the  two  dimensional  symmetry  of  flames  and  to  improve  uniformity 
of  the  fuel  flow.  In  the  original  two-hole  case,  the  fuel  flow  was  under  a  significant  end  boundary 
condition  effect  from  the  post  combustion  stream  since  only  one  fuel  hole  in  each  side  of  the  bluff 
body  existed.  Two-dimensional  symmetry  from  the  configuration,  therefore,  was  not  achievable. 
When  the  number  of  holes  was  increased  to  ten  (five  on  each  side),  it  was  found  that  the 
flame  uniformity  from  each  hole  was  problematic  due  to  a  pressure  gradient  along  the  1/4” 
fuel  delivery  pipe;  this  was  alleviated  significantly  if  we  used  only  six  holes.  Thus,  the  final  6 
hole  configuration  obtained  by  deleting  4  holes  at  the  boundaries  from  the  10  hole  configuration 
satisfactorily  provided  both  the  two  dimensional  symmetry  in  the  central  region  of  the  bluff 
body  and  good  flame  uniformity.  All  experiments  described  here,  therefore,  have  been  carried 
out  using  the  6  hole  configuration. 

4.1.2  Implementation  of  Water  Cooled  Grid  and  Turbulence  Grid 

The  overall  experimental  setup  with  a  newly  added  (from  the  setup  of  Phase  I)  water-cooled  grid 
and  a  turbulence  grid  is  provided  in  figure  1.  While  the  main  burner  is  unchanged  from  Phase 
I,  the  water-cooled  copper  grid  which  has  a  similar  geometry  to  the  condenser  of  a  refrigerator 
is  placed  on  the  top  of  the  main  burner.  The  use  of  the  water  cooled  grid  eliminates  the 
need  for  additional  air /nitrogen  injection  into  the  post  combustion  gas  (as  used  throughout 
Phase  I),  which  added  to  the  non-uniformity  of  the  velocity  and  concentration  fields  of  the  post 
combustion  stream.  The  turbulence  grid  is  placed  ~  1"  downstream  of  the  top  surface  of  the 
water-cooled  grid  (see  Fig.  1  for  detailed  dimensions).  Two  opening  spacings  of  the  grid  (25% 
and  36%  opening  areas)  and  two  distances  between  the  turbulence  grid  and  the  bluff  body,  0.5" 
and  2",  have  been  tested  to  optimize  the  level  of  turbulence.  The  detailed  geometries  of  the 
water  cooled  grid  and  the  turbulence  grid  are  shown  in  Figs.  2(a)  and  2(b),  respectively.  Water 
at  room  temperature  flows  through  a  1/4"  (OD),  serpentine-shaped  copper  tube  which  has  11 
turns  within  a  15"  by  18"  rectangular  region.  Numerous  1/16"  (OD)  copper  wires  are  welded  on 
the  copper  tube  to  improve  heat  conduction  along  the  streamwise  direction  of  water  flow.  The 
resulting  temperature  profile  along  the  tube,  therefore,  can  be  more  uniform.  The  turbulence 
grid  (see  Fig.  2(b))  is  composed  of  two,  top  and  bottom,  layers  of  1/2"  OD,  8"  long  ceramic 
(99.9%  A1203)  tube  arrays.  Each  array  is  formed  by  6,  7,  or  8  parallel-placed  ceramic  tubes, 
which  depends  on  the  size  of  opening  areas  (7  for  top  layer  of  25  %  opening;  8  for  bottom  layer 
of  25  %  opening;  and  6  for  top  and  bottom  layers  of  36  %  opening).  The  top  and  bottom  arrays 
are  placed  perpendicular  to  each  other  such  that  the  shape  of  the  open  area  becomes  square. 
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D1:  distance  from  tfie  burner  top  surface  to  the  center  of  water  grid  (1  1 


02:  distance  from  the  center  of  water  grid  to  the  center  of  bottom  ceramic  pipe  array  (1  1/8") 
D3:  distance  from  the  center  of  water  grid  to  the  center  of  top  ceramic  pipe  array  (1  5:8”) 


Figure  1:  Schematic  of  overall  experimental  setup 


The  distances  between  the  centers  of  adjacent  tubes  are  0"  (25  %  opening)  and  1.25"  (36  % 
opening),  respectively. 

Figure  3  shows  sample  PIV  results  to  compare  the  turbulent  kinetic  energyT  \  {uf2+ 1/2},  where 
u  is  the  horizontal  and  v  is  the  vertical  velocity  component ,  respectively,  for  three  different  cases: 
Without  the  turbulence  grid  (Fig.  3a),  with  the  turbulence  grids  with  25  %  opening  area  (Fig. 
3b)  and  36  %  opening  area  (Fig.  3c).  For  all  cases  shown  in  Fig.  3,  50  individual  images 
are  ensemble- aver  aged  and  the  water  cooled  grid  is  always  implemented.  Color  coding  of  each 
figure  represents  the  intensity  of  the  turbulent  kinetic  energy  (high  turbulent  kinetic  energy 
in  dark  red  regions  and  low  in  dark  blue  regions).  In  all  three  cases,  initial  flow  conditions 
(MX  65  equivalence  ratio  in  the  inlet  of  the  main  burner,  ^3.5  m/s  free  stream  speed,  -^1200 
K  free  stream  temperature,  no  bluff  body)  and  imaging  region  (~70  mm  x  70  mm  square 
w'hose  bottom  line  is  located  at  5"  above  the  top  surface  of  the  turbulence  grid)  are  identical. 
As  shown,  the  no-grid  case  has  average  turbulent  kinetic  energy  of  fV  0.1  (m/s)2  (or  Unarm  « 
(• u 2  4-  v2)0-5 /vfreest ream  —  13%)  while  25%  and  36%  opening  cases  have  ~0.17  m/s2  and  ~0.22 
m/s2,  respectively  (or  Un^m  —^17  %  and  ^19  %,  respectively).  Therefore,  one  can  conclude 
that  the  36  %  opening  area  grid  has  a  higher  capability  in  generating  more  intense  turbulence 
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Figure  2:  Detailed  dimensions  of  (a)  water  cooled  and  (b)  turbulence  grids 
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Figure  3:  Turbulence  intensities  in  the  presence  of  (a)  no  grid,  (b)  25  %  opening  area  grid  and 
(c)  36  %  opening  area  grid 


than  the  25  %  opening  area  grid. 

Other  sample  PIV  results  obtained  in  the  presence  of  the  bluff  body  are  shown  in  Fig.  4. 
Figures  4a  and  4b  show  the  results  in  the  absence  and  presence  of  the  turbulence  grid  (36  % 
opening  area),  respectively.  Initial  Sow  conditions  axe  identical  wdth  those  of  Fig.  3,  but  the  PIV 
imaging  region  is  adjusted  to  1.5"  higher  than  that  of  Fig.  3  to  investigate  the  possible  decay  of 
turbulence  intensities  further  downstream  than  the  case  of  Fig.  3.  While  the  wake  region  due 
to  the  presence  of  the  bluff  body  is  well  captured  by  the  dark  red  region  in  both  figures,  one  can 
observe  that  the  case  with  the  turbulence  grid  (Fig.  4b)  still  has  higher  freestream  turbulence 
intensities  {unorm  =  ^23  %)  than  that  without  the  grid  ( unorm  =  ^19  %), 

4*1*3  Effect  of  OH  Radical  Addition  Though  Hydrogen  Peroxide  (H2Q2)  Injection 

The  injection  of  H2Q2  into  the  post  combustion  gas  is  carried  out  with  the  setup  shown  in 
Fig.  5.  The  steady  flow  of  liquid  H2O2/H2O  solution  (pure  water,  30  vol.%  or  50  vol.%  H2O2 
solution  in  w^ater)  is  implemented  by  a  step  motor/syringe  assembly  placed  ^2  m  upstream  of 
the  main  burner,  while  the  liquid  flows  through  four  1/8"  OD,  2  m  long  teflon  tubes  which 
connect  the  syringe  and  four  central  tubes  in  the  top  array  of  the  turbulence  grid.  The  total 
flowrate  of  the  H2O2/H2O  solution  provided  by  the  step  motor/syringe  assembly  is  ^20  seem 
(liquid)  for  each  experiment.  The  solution  is  in  liquid  phase  until  it  reaches  the  inlets  of  the 
ceramic  tubes,  but  is  rapidly  vaporized  when  it  flows  through  the  inside  of  the  ceramic  tubes. 
The  resulting  H2O2/H2O  vapor  is  injected  into  the  post  combustion  gas  through  five  —  1/8" 
x  1/16"  rectangular  shaped  slits  in  each  of  the  four  ceramic  tubes.  The  inter-spacing  between 
adjacent  holes  on  one  ceramic  tube  is  1".  Flame  stabilities  with  the  H2O2/H2O  injection  are 
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Figure  4:  Turbulence  intensities  in  the  presence  of  bluff  body,  (a)  no  grid,  (b)  36  %  opening 
area  grid.  Imaging  region  is  1,5”  higher  than  that  of  figure  3 


quantified  by  investigating  the  liftoff  height  of  an  afterburning  methane  flame.  As  performed 
in  Phase  I,  CH  chemiluminesence  of  the  fiamebase  is  detected  by  a  100  s  gate  width  intensified 
CCD  camera  (Princeton  Instrument,  PPM  AX)  to  determine  instantaneous  fiamebase  locations. 
The  liftoff  height  is  determined  by  an  average  over  200  ensembles  of  the  instantaneous  fiamebase 
location. 

Figure  6  shows  a  sample  result  of  the  flame  liftoff  height  measurements  under  four  different 
amounts  of  the  H2O2/H2O  injections:  no  injection,  100  %  water  (20  seem  total),  30  %  H2Q2  in 
70  %  water  (20  seem  total),  and  50  %  H2O2  in  50  %  water  (20  seem  total).  We  are  limited  to 
50%  H2O2  due  to  safety  considerations.  Flow  conditions  axe  identical  with  those  of  Fig,  4  except 
for  the  injections.  As  shown,  the  flame  stability  decreases  with  increasing  water  concentration 
of  the  injection.  In  comparison  with  the  nounjection  case,  for  example,  100  %  water  injection 
decreases  the  flame  stability  by  ^  25  %.  However,  the  increasing  H2O2  portion  in  the  injection 
reverses  the  propensity  such  that  the  flame  stability  of  the  50  %  H2G2  addition  is  comparable 
with  that  of  the  no- inject  ion  case.  It  is  noteworthy  that,  based  on  our  preliminary  PREMIX 
calculation  (not  shown  here),  more  than  99.99  %  of  H2O2  is  dissociated  to  form  H2O  at  this 
temperature  (~  1200  K)  within  less  than  1  ms  while  ^  IQ-4  mole  fraction  OH  radicals  still 
linger.  Therefore,  increasing  flame  stability  with  increasing  amounts  of  H2O2  addition  shown 
in  this  result  provides  plausible  evidence  on  the  central  role  of  the  OH  radicals  on  the  flame 
stability. 

In  summary,  the  results  shown  here  provide  a  configuration  for  experimental  studies  to 
generate  higher  turbulence  intensity  and  higher  OH  radical  concentration  in  the  post  combustion 
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Figure  5:  Turbulence  intensities  in  the  presence  of  bluff  body,  (a)  no  grid,  (b)  36  %  opening 
area  grid.  Imaging  region  is  1.5”  higher  than  that  of  figure  3 
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Diluted  species 

Figure  6:  Flame  liftoff  heights  under  various  amounts  of  H2O2/H2O  injections 


stream.  Future  studies  can  allow  a  broader  range  of  the  turbulence  intensity  and  the  OH  radical 
concentration,  detailed  measurements  of  the  bluff-body  shedding  frequencies  and  spectra  and 
simultaneous  and  instantaneous  CH  PLIF/PIV  measurements  using  ethylene  fuel  for  greater  CH 
signals.  Such  results  will  provide  a  more  comprehensive  dataset  which  will  allow  a  more  careful 
determination  of  the  effect  of  radical  composition  and  the  underlying  flow  field  on  bluff-body 
flame  stability. 

4.2  Modeling  and  Computation 

4.2.1  Validation  and  Damkohler  Number  Similarity 

Flame-holders  provide  two  flame  stability  mechanisms.  As  the  flame-holder  modifies  the  flow 
around  it,  a  wake  is  generated  aft  of  the  flame-holder  base.  This  wake  provides  sufficient  res¬ 
idence  time  and  chemical  heat  release  needed  to  ignite  the  fresh  gas.  In  addition,  the  shear 
layer  generated  around  the  wake  provides  a  turbulent  mechanism  to  mix  the  fresh  air  with  the 
combustion  products  from  the  recirculation  zone  [1,  2].  Finally,  as  the  flame  is  stabilized  past 
the  flame-holder  base,  combustion  products  at  the  end  of  the  recirculation  zone  are  entrained 
into  the  wake  to  provide  a  self-ignition  mechanism. 

The  inflow  mixture  to  the  augmentor  usually  has  high  speed  and  high  temperature.  The 
typical  inflow  temperature  and  velocity  for  augmentors  range  between  700— 1000  K  and  100—200 
m/s  respectively  [3].  Under  these  conditions,  the  experimental  assessment  of  such  configurations 
is  complex  and  expensive,  and  it  becomes  practically  challenging  to  investigate  the  flow  and 
combustion  features.  The  current  work  investigates  the  possibility  of  utilizing  similarity  to  test 
such  configurations.  Here,  we  will  assess  the  Damkohler  number  similarity,  which  implies  that 
in  the  high  Reynolds  number  limit,  the  flow  and  the  flame  characteristics  in  a  high-speed/fast 
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chemistry  case  should  be  similar  to  a  low-speed/slow  chemistry  case,  as  long  as  the  Damkohler 
number  and  all  the  momentum  ratios  is  the  same. 

Here,  three  numerical  simulations  are  performed  using  large-eddy  simulations.  The  chemical 
source  term  is  modeled  using  the  hamelet  progress  variable  (FPV)  approach.  The  configuration 
is  for  a  bluff-body  stabilized  Same  in  a  vitiated  cross- Sow.  Methane  fuel  is  injected  in  the  cross- 
stream  direction  from  the  bluff-body  surface.  The  geometry  of  the  bluff  body  is  represented  in 
the  simulation  by  an  immersed  boundary  technique  on  a  structured  grid.  The  combination  of 
immersed  boundary  with  structured  codes  gives  the  advantage  of  simulating  complex  geometries 
with  lower  computational  cost  and  high  accuracy  because  of  the  staggered  variable  arrangement. 
For  constant  Reynolds  number,  the  inSow  temperature  and  chemical  composition  are  varied  to 
change  the  mutual  effect  between  mixing  and  chemical  time  scales,  and  hence  the  Damkohler 
number.  The  inflow  chemical  variation  and  the  Damkohler  number  are  found  to  affect  both  the 
flame  stabilization  and  the  lift-off  characteristics.  The  stagnation  point  position  in  the  wake  of 
the  bluff  body  and  the  lift-off  height  are  monitored  for  both  test  cases  to  assess  the  Damkohler 
number  similarity.  This  section  is  organized  as  follows:  We  will  first  present  the  numerical 
and  modeling  approaches  used,  the  non-dimensional  form  of  the  governing  equations,  and  the 
definition  of  the  characteristic  Damkohler  number.  Next,  we  will  show  the  numerical  setup  and 
the  geometry  details.  Finally,  the  results  will  be  discussed. 

Numerical  Formulation  -  All  simulations  reported  here  were  performed  with  the  Stanford 
NGA  code.  The  numerical  scheme  used  here  is  based  on  a  high  order  conservative  finite  dif¬ 
ferences  [4,  5j.  The  governing  equations  are  in  the  low  Mach  number  formulation.  The  scheme 
is  staggered  in  space  and  time  [6].  The  staggering  allows  lor  secondary  conservation.  Hence, 
no  artificial  damping  is  required  and  the  scheme  can  allow  for  higher  order  accuracy  [5],  How¬ 
ever,  in  the  following  section  we  will  focus  only  on  the  details  of  the  combustion  model  and  the 
Damkohler  number  similarity  analysis. 

By  applying  an  implicit  top  hat  filter  denoted  by  an  overline  to  the  Navier  Stakes  equations 
and  using  the  Favre  decomposition,  defined  for  a  general  variable  <b  as  $  =  fxp/p,  the  continuity 
and  the  momentum  equations  can  be  written  as: 

dp  dpnn 
di  dxt 

dm  dpUjUj  _  .  dp  _  dr ij_ 

&t  dxj  dxi  dxj  dxj 

Here,  Ui  is  the  filtered  velocity,  p  is  the  filtered  density,  P  is  the  filtered  pressure,  and  the 
filtered  shear  stress  is  written  +  The  summation  over  repeated 

indices  is  implied  and  fi  is  the  dynamic  viscosity.  In  the  above  equations,  the  influence  of  the 
statistical  effect  of  the  subgrid  scales  on  the  resolved  scales  is  represented  by  the  sub-filter  stresses 
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Tij  —  p  (UiUj  -UjUj).  The  sub-filter  stresses  are  modeled  with  an  eddy  viscosity  assumption  [7): 


Tij  —  “2mT  ~  , 


(2) 


where  Stj  =  \  is  the  velocity  strain  rate,  and  q2  is  the  trace  of  the  sub-filter  stress 

tensor.  The  turbulent  eddy  viscosity  px  and  q2  are  modeled  as: 


W=Crf&2\S\ 

q2  =  2CkpA2\S\2.  (3) 

In  these  equations  |5|=  ^2 IhjSij  and  and  Cjt  are  computed  dynamically  based  on  the 
approach  introduced  by  Moin  et  ah  [7]*  The  density  p  is  calculated  from  the  state  equation 
described  in  the  next  section.  Under  the  low  Mach  number  approximation  [8],  the  acoustic 
modes  are  decoupled  from  the  entropy  and  the  vorticity  modes.  This  approximation  allows  for 
a  larger  and  more  stable  integration  time  step  for  Eqs,  (1),  A  semi-implicit  iterative  technique 
is  employed  to  integrate  the  above  system  of  equations  using  the  fractional  step  method*  More 
details  about  the  above  numerical  approach  can  be  found  elsewhere  [6,  5]  and  hence  are  avoided 
here  for  brevity*  The  following  section  will  describe  the  employed  combustion  model. 

Flamelet  Progress  Variable  Approach  -  The  fiamelet /progress  variable  (FPV)  approach 
developed  by  Pierce  and  Moin  [9]  is  based  on  the  flamelet  concept  [10].  This  approach  assumes 
that  the  chemical  time  scales  are  fast  enough  such  that  the  reactions  occur  in  thin  layers  around 
the  stoichiometric  conditions  [11].  The  flamelet  approach  relates  the  species  mass  fractions  and 
the  energy  to  the  mixture  fraction  through  the  fiamelet  equations  [12]: 


0$  XdH  . 

P  dt 


(4) 


where  Z  is  the  mixture  fraction,  $  =  (T,yi)r  is  a  vector  that  contains  the  species  mass  fractions 
yi  and  the  temperature  T .  u)  is  the  source  term  and  x  is  the  scalar  dissipation  rate,  which 
represents  the  inverse  of  the  diffusion  characteristic  time  scale  or  equivalently  the  fluid  influence 
on  the  nan-conserved  scalars.  In  Eq.  (4)  the  scalar  dissipation  rate  is  given  by  [13] 


X  =  2D\VZ\2.  (5) 

In  the  FPV  model,  a  new  reactive  scalar  A  is  introduced.  As  a  result,  the  state  vector  is 
given  as  (Z,  A).  Accordingly,  this  formulation  assumes  that  the  reactive  scalars  are  a 

function  of  Z  and  A*  Then 

noo 

(Z,  A)  P  (Z,  A;  x,  £)  dXdZ .  (6) 

_ 

The  A-parameter  is  represented  by  the  progress  variable  C  that  can  be  a  combination  of  certain 
species  products,  temperature,  enthalpy  or  any  scalar  that  can  represent  the  reactive  state 
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uniquely.  The  progress  variable  here  is  taken  as  the  sum  of  the  mass  fractions  of  CO2,  CO, 
H2O,  and  Bfe.  Then,  also  C  can  be  obtained  from  Eq.  6,  and  A  can  be  eliminated  from  the 
problem.  The  passive  scalar  marginal  PDF  P  ( Z )  is  assumed  to  be  a  beta  distribution  [9] ,  and 
the  reactive  scalar  PDF  as  delta  function.  Since  the  PDF  distribution  is  characterized  by  the 
scalar  mean  and  the  variance,  the  flamelet  tabulation  will  take  this  form: 

$  =  ^(z,F2,c).  (7) 

Therefore,  three  parameters  have  to  be  computed,  namely  the  filtered  mixture  fraction  Z,  the 
sub-filter  scalar  variance  Z "2,  and  the  filtered  progress  variable  C.  Here  the  scalar  variance 
Z"2  is  evaluated  using  an  algebraic  model  [14].  The  model  assumes  homogeneity  and  local 
equilibrium  for  the  sub-grid  scales  and  is  given  as 


pZ"2  =  Cztfp\VZ\2,  (8) 

where  the  mixture  fraction  coefficient  Cz  is  computed  dynamically.  The  second  parameter  is 
the  filtered  mixture  fraction  Z,  for  which  a  transport  equation  is  solved  as 

^  +  V-(pu2)  -V-  (p{Dz  +  DT)VZ)  .  (9) 

The  mixture  fraction  diffusivity  D%  is  obtained  from  Eq.  (6),  and  Dj  is  the  turbulent  eddy 
diffusivity  given  by  [14] 


Dt  =  CtA?\S\,  (10) 

where  Ct  is  computed  in  a  dynamic  procedure  [7].  Finally,  the  filtered  progress  variable  is 
obtained  by  solving  the  transport  equation  of  C  as: 

+  V  •  (puC)  =  V  •  [p{Dc  +  DT)  VC)  +  Wc ,  (11) 

where  Dq  is  the  progress  variable  diffusivity,  and  Wc  =  p{mc/p)  is  the  Favre  averaged  net  mass 
production  rate  of  C,  me  is  the  unfiltered  net  production  rate  in  (kg/m3/sec).  In  the  above 
formulation,  given  the  mean  and  the  variance,  the  subgrid  fluctuations  effect  on  the  filtered  state 
vector  $  is  represented  by  the  integration  over  the  joint  PDF  of  Z  and  C  expressed  in  Eq.  (6) 
[15].  Next  we  will  reformulate  the  governing  equations  in  the  non-dimensional  form  to  show  the 
dependence  on  the  Damkohler  number  and  the  similarity  conditions. 

Non-Dimensional  Governing  Equations  -  To  derive  the  non-dimensional  governing  equa¬ 
tions,  the  following  relations  are  defined,  where  the  subscript  (00)  refers  to  the  reference  values 
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at  the  vitiated  inflow  conditions  and  the  superscript  (*)  refers  to  the  non-dimensional  values: 
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Here  the  reference  length  scale  is  chosen  as  the  bluff-body  diameter  Db.  In  the  following  analysis, 
the  inflow  progress  variable  diffusivity  Dc,oo  and  the  inflow  mixture  fraction  diffusivity  Dz,<x> 
are  fixed  to  a  single  value  that  corresponds  to  the  maximum  inflow  diffusivity.  Another 

important  parameter  r  =  \j  is  the  momentum  ratio  between  the  fuel  jet  (subscript  j)  and 
the  vitiated  oxidizer  stream  (subscript  oo).  The  scalar  dissipation  rate  x  is  non-dimensionalized 
by  its  value  at  the  quenching  conditions  Xq-  Note  that  the  mixture  fraction  as  a  normalized 
quantity  needs  no  further  normalization. 

By  substitution  of  Eqs.  (12)  in  Eqs.  (1)  (and  dropping  the  *),  we  get  the  following  non- 
dimensional  form: 


dp  dpuj 
9t  dxi  * 

dpuj  dpujUj  8P  ,  1  (d&g\  dry 

dt  dXj  dXi  dXj  )  dxj  '  (  ’ 

where  the  vitiated  inflow  Reynolds  number  is  given  by  Re^  =  In  Eq.  (13),  the  subgrid 

models  presented  by  Eq.  (3)  are  used  to  non-dimensionalize  the  sub-filter  stresses  ry.  Following 
the  same  procedure  with  Eqs.  (9)  and  (11),  we  get  the  filtered  progress  variable  and  the  filtered 
mixture  fraction  equations  in  the  non-dimensional  form  as 


where  Scoo  =  is  the  vitiated  inflow  Schmidt  number  for  the  scalar  progress  variable  and 

the  mixture  fraction  and  the  non-dimensional  source  term  and  the  Damkohler  number  Da^  will 
be  introduced  below.  Further  treatment  for  the  progress  variable  source  term  will  be  discussed 
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in  the  next  subsection.  Finally,  the  state  relationship  for  the  thermodynamic  and  chemical 
variables  is  computed  from  the  non-dimensional  equation  of  state: 

=  .  (16) 

For  the  same  non-dimensional  quantities  (i.e,  Re^,  Sc^,  and  Dac),  Eqs.  (13)-(16)  will  have  the 
same  solution  for  similar  boundary  conditions.  The  inflow  boundary  conditions  for  the  cross 
flow  and  the  jet  are  uniform  with  bulk  velocity  UM  and  Uj  respectively.  Therefore,  the  inflow 
boundary  conditions  for  the  vitiated  cross  flow  in  the  non-dimensional  form  will  be: 


utoo  =  L0 

u2,oo  =  0-0 
u3,oo  =  0-0 

c~  =  1.0 

Z  =  0.0. 

Similarly,  the  non-dimensional  jet  inflow  conditions  will  be: 


u*,j  =  0-0 
u2,j  =  0-0 


C*  =  0.0 
Z  =  1.0. 


(17) 


(18) 


where  r  =  \i is  the  momentum  ratio,  and  p*  =  •&-.  The  convective  condition  is  used  for 

V  POQ^OO  J  poc 

the  outflow  boundary.  This  condition  can  be  written  for  a  general  variable  ip  as: 


dip  dip 
9t  +  a  dxn 


=  0, 


(19) 


where  a  is  the  convection  speed  at  the  outflow  cross-section,  and  xn  is  the  coordinate  in  the 
direction  of  the  outward  normal  at  the  boundary.  In  the  non-dimensional  form  this  equation 
can  be  written  as: 


dZi  *  dui 
dt  '  dx„ 

dC*  ,dC* 
dt  a  dxn 
dZ  ^  *dZ_ 
dt  '  a  dxn 


=  0 
■  =0 
=  0. 


(20) 


The  above  equations  show  that  given  two  different  test  cases,  with  the  same  momentum  ratio 
(r)  and  the  same  non-dimensional  reference  quantities  (i.e.  ReM,  Scoo,  and  Dac),  the  set  of 
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Eqs.  (13)-(15)  will  have  similar  boundary  conditions  and  similar  solution*  However ,  in  the  next 
subsection  the  source  term  in  Eq.  (15)  will  be  discussed  in  more  details. 

Damkohler  Number  Similarity  -  The  characteristic  Damkohler  number  is  defined  as  the 
ratio  between  the  characteristic  flow  or  mixing  time  scale  to  the  characteristic  chemical  time 
scale 

Dac  =  —  (21) 

^chem 

where  the  flow  time  scale  is  defined  as  Tflow  =  based  on  the  bluff-body  width  D &  =  25.4  mm 
and  the  cross  flow  velocity  U^.  The  chemical  time  scale  can  be  defined  by  many  different  meth¬ 
ods.  For  single  step  chemistry,  it  is  typically  defined  as  the  chemical  source  term  divided  by  a 
reactant  concentration*  In  a  previous  work.  Pitsch  and  Fedotov  [16]  derived  an  expression  for  the 
non-dimensional  source  term  for  a  single  step  chemistry  They  showed  that  the  source  term  and 
equivalently  the  chemical  time  scale  are  a  function  of  unique  non-dimensional  quantities  that 
depend  on  the  frequency  factor  of  the  one-step  mechanism  and  the  chemical  and  thermal  bound¬ 
ary  conditions.  However,  detailed  chemical  mechanisms  describe  many  combustion  phenomena, 
such  as  auto-ignition,  premixed  flame  propagation,  and  non-premixed  combustion.  Because  of 
the  different  interactions  of  these  processes  with  molecular  transport,  different  chemical  reac¬ 
tions  might  be  important  in  these  processes,  and  there  is  no  unique  definition  of  the  relevant 
chemical  time  scale.  Each  chemical  reaction  has  its  own  time-scale  and  Damkohler  number.  If  a 
similar  solution  is  to  be  achieved*  for  instance,  at  a  higher  flow  velocity,  which  implies  a  smaller 
flow  time  scale,  also  the  chemical  time  scales  need  to  be  reduced  to  achieve  the  same  Damkohler 
numbers*  Typical  ways  to  change  the  chemical  time  scales  are  changes  in  the  temperature,  the 
equivalence  ratio,  or  dilution*  It  is  obvious  that  not  all  chemical  time  scales  can  be  changed  in 
the  same  way.  Therefore,  the  choice  of  the  adjusted  chemical  time  scale  should  be  related  to 
the  physics  of  the  dominating  combustion  phenomenon.  For  example,  in  premixed  flames,  the 
flame  response  based  on  the  flame  thickness  and  flame  propagation  speed  is  more  relevant.  For 
non-premixed  flames,  the  chemical  time  scale  should  represent  the  mutual  interaction  between 
chemistry  and  mixing*  The  following  analysis  shows  that  for  the  problem  in  consideration  the 
vitiated  inflow  temperature,  for  a  given  inflow  speed,  has  a  strong  effect  on  the  available  range 
of  chemical  time  scales  and  on  the  ignition  and  quenching  characteristics*  As  shown  in  Fig.  7,  as 
the  vitiated  inflow  temperature  increases,  the  production  rate  of  certain  species  is  accelerated, 
and  the  mixture  is  ignited  at  a  faster  chemical  time  scale*  This  increases  the  strain  rate  and  the 
mixture  fraction  gradient,  and  allows  for  a  higher  quenching  scalar  dissipation  rate  as  well.  As 
the  inflow  temperature  increases  more,  around  =  1700  K  the  mixture  is  self-ignited  and  the 
lower  unstable  branch  disappears. 

Figure  8  shows  the  non-dimensional  source  term  versus  the  non-dimensional  dissipation  rate 
for  different  vitiated  inflow  temperatures.  The  source  term  in  Eq.  (15)  will  be  non-dimensionlized 
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Figure  7:  The  effect  of  the  vitiated  hot  flow  temperature  (Tjn/)on  quenching  and  ignition  char¬ 
acteristics  of  the  flamelets  at  the  stoichiometric  conditions 


with  the  maximum  progress  variable  production  rate  Wc.moi  as 


WA  = 


WC 

Wc^rnax 


The  corresponding  characteristic  chemical  time  scale  is  then 

_  _  Poo  Coo 

'chem  rir  i 

Wc,max 

and  the  Damkohler  number  becomes 

Da^  —  C.max 

UooPoo  Cqo 


(22) 


(23) 


(24) 


It  is  obvious  that  the  minimum  chemical  time  scale  appears  approximately  at  stoichiometric 
mixture  fraction.  To  illustrate  this  observation,  Fig.  9  shows  a  scatter  plot  of  the  inverse  of 
the  local  chemical  time  scale  with  the  mixture  fraction.  The  minimum  chemical  time  scale  that 
corresponds  to  the  highest  production  rate  is  located  at  the  stoichiometric  condition.  As  will  be 
shown  next,  the  value  of  the  production  rate  at  quenching  is  slightly  lower  than  the  maximum 
value. 

Further  insight  can  be  achieved  by  plotting  the  three  dimensional  S-shaped  curve  in  Fig.  10, 
which  shows  the  maximum  temperature  and  the  chemical  time  scale  as  function  of  the  dissipation 
rate.  The  curve  shows  that  the  upper  and  lower  branch  have  different  time  scales,  which  is  even 
more  clear  in  Fig.  10b.  For  the  burning  branch  the  chemical  time  scale  is  smaller.  As  also  shown 
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Figure  8:  The  variation  of  the  non-dimensional  progress  variable  production  rate  with  the  non- 
dimensional  scalar  dissipation  rate 


Figure  9:  Scatter  plot  for  the  inverse  of  the  local  chemical  time  scale  with  the  mixture  fraction 
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(b)  Projection  in  the  x-  time  inverse 
plane 


Figure  10:  Three  dimensional  S-shaped  curve 


in  Fig.  10b,  the  value  of  the  chemical  time  scale  at  the  quenching  point  is  slightly  higher  than 
the  minimum  value. 

It  should  be  noted  again  that  here  the  characteristic  chemical  time  scale  is  defined  as  the 
lowest  possible  time  scale  at  stoichiometric  conditions.  This  time  scale  corresponds  to  the 
maximum  Damkohler  number  for  a  given  flow  condition. 

Numerical  Setup  -  The  bluff-body  geometry  and  the  grid  after  refinement  are  shown  in 
Fig.  11.  Methane  jets  at  300  K  are  injected  from  two  circular  ports  0.5  inch  upstream  the  base 
of  the  bluff  body.  A  flow  of  vitiated  hot  mixture  is  injected  in  the  stream-wise  direction.  The  jet 
and  the  vitiated  mixture  mass  flow  rates  are  varied  to  change  the  effective  flow  time  and  hence 
the  Damkohler  number.  The  bluff  body  is  4.5  inch  long  with  a  1  inch  width  rectangular  body 
and  a  1/2  inch  radius  rounded  nose. 

Since  the  current  work  is  devoted  to  the  validation  and  assessment  of  the  similarity  solution, 
three  different  test  cases  are  chosen  here.  The  flow  conditions  for  the  three  test  cases  are 
summarized  in  Table  1. 

The  conditions  of  the  first  test  case  match  the  setup  of  the  Stanford  vitiated  flow  experi¬ 
ment,  while  in  the  other  two  cases  the  flow  and  the  chemical  conditions  are  varied  such  that  the 
Damkohler  number  is  kept  the  same.  The  first  test  case  represents  low-speed/slow  chemistry, 
while  the  second  and  the  third  case  represent  high-speed/fast  chemistry.  The  inflow  composi¬ 
tion  is  computed  as  the  equilibrium  composition  of  methane/air  at  the  inflow  equivalence  ratio 
and  temperature.  These  values  are  then  used  as  an  input  to  the  flamelet  equations  boundary 
conditions. 

Results  -  The  instantaneous  recirculation  zones  (RZ)  and  flame  structures  for  all  three  cases 
are  shown  in  Fig.  12.  The  three  flames  show  similar  flame  and  RZ  structure.  The  recirculation 
zone  is  found  to  extend  over  the  bluff-body  base  and  is  surrounded  by  the  flame  surface.  The 
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Figure  11:  Computational  mesh  for  LES  of  the  flame  holder  geometry 


Case# 

Uoo 

T 

Tootif) 

Moo 

Dao 

*00 

Tfloxu 

T. chetn 

Dac 

Reoo 

Scqc 

Case  1 

3.6 

2.23 

1.3 

1100 

1.6E-5 

O.lE-3 

0.67 

0.0071 

0.0011 

6.429 

5715 

0.16 

Case  2 

7.8 

4.43 

1.3 

1300 

3.45E-5 

0.22E-3 

0.67 

0.00327 

0.00051 

6.429 

5715 

0.16 

Case  3 

22.7 

9.75 

1,3 

1700 

l.E-4 

4.23E-3 

0.67 

0.00118 

1.74E-4 

6.429 

5715 

0.16 

Table  1:  Flow  test  conditions,  where  <3>oo  and  Too  are  the  equivalence  ratio  and  the  temperature 
of  the  vitiated  hot  stream  respectively,  fi0 0  is  the  inflow  kinematic  viscosity  in  (m2/s)?  D ^  is 
the  inflow  molecular  mass  diffusivity  (m2/s),  Uoc  is  the  vitiated  air  bulk  velocity  in  (m/s),  Uj 
is  the  fuel  jet  velocity  in  (m/s),  7fiow  is  the  characteristic  flow  time  in  seconds,  rChem  is  the 
characteristic  chemical  time  in  seconds,  Dac  is  the  characteristic  Damkohler  number,  and  Reoo? 
and  Scoo  are  the  vitiated  hot  flow  Reynolds  number  and  Schmidt  number  respectively. 
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(a)  Coae(l)  (b)  Case(2)  (c)  Case(3) 


Figure  12:  The  iso-surface  of  negative  axial  velocities  (grey)  combined  with  the  temperature 
iso-surface  colored  by  the  progress  variable 

effect  of  the  outer  shear  layer  is  shown  to  wrinkle  the  flame  surface,  causing  shedding  of  mixed 
burned/unburned  gases  towards  the  outflow.  More  discussion  will  be  provided  below. 

Figure  13  shows  the  mean  non-dimensional  axial  velocity  distributions  at  different  heights 
above  the  bluff-body  base.  The  RZ  extends  to  a  distance  of  around  50  mm  above  the  base. 
The  wake  extends  about  20  mm  in  both  directions  of  the  centerline.  Experimental  data  are 
available  at  three  locations  and  are  compared  with  the  simulation  results.  Good  agreement  is 
noticed  between  the  LES  and  the  experimental  data.  However,  the  LES  data  show  a  symmetrical 
solution,  while  the  experimental  data  do  not.  This  asymmetry  in  the  experimental  results  is 
due  to  the  fact  that  the  inflow  field  is  not  perfectly  uniform.  Most  importantly  here,  the  non- 
dimensional  curves  collapse  onto  each  other,  which  confirms  the  flow  similarity  of  the  three 
flames. 

A  comparison  of  the  normalized  mean  temperature  is  shown  in  Fig.  14.  The  three  test  cases 
exhibit  the  same  flame  structure.  However,  at  the  first  location,  the  case  (3)  flame  shows  a 
small  deviation  at  the  edge  of  the  shear  layer.  At  this  location,  the  effect  of  the  unsteadiness 
due  to  the  RZ  oscillation  contributes  to  this  deflection  as  will  be  discussed  later.  The  above 
results  show  that  the  mean  flame  and  RZ  structure  exhibit  a  similarity  solution  under  the  given 
conditions  for  the  three  studied  flames. 

As  the  studied  flames  show  similar  RZ  structure,  they  also  show  similar  flame  stability 
features.  The  flame  stabilization  point  is  defined  as  the  first  point  where  the  flame  and  the  flow 
speed  match.  Since  this  flame  is  wake  stabilized,  the  wake  structure  controls  the  flame  stability 
characteristics.  Furthermore,  as  long  as  the  chemical  time  scale  is  fast  enough  compared  with 
the  flow  time  scales,  the  flame  will  stabilize  and  will  not  blowoff.  The  faster  the  chemical  time 
scale,  the  higher  the  flow  speed  needed  to  stabilize  the  flame  at  the  same  height.  This  finding 
ascertains  the  existence  of  Damkohler  number  similarity  between  the  three  flames.  To  further 
asses  the  validity  of  the  current  approach,  the  flame  height  above  the  base  is  measured  and 
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Figure  13:  Mean  non-dimensional  axial  velocity  ( jj2-)  at  different  heights  X  above  the  bluff-body 
base;  *  experimental  PIV  data,  (— )  case  (1)?  (*)  case  (2)*  (o)  case  (3) 
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Figure  14:  Mean  normalized  temperature  (0  =  )  at  different  heights  above  the  bluff- 

body  base,  (— )  case  (1),  {*)  case  (2),  (o)  case  (3).  T$t  is  the  maximum  burning  temperature  at 
stoichiometric  conditions 
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compared  with  experimental  data.  By  analogy  with  the  experiment,  the  flame  height  is  defined 
as  the  first  point  above  the  base  (at  the  median  plane) ,  where  the  stoichiometric  conditions  are 
found.  The  flame  height  is  computed  as  the  average  over  the  entire  instantaneous  snapshots. 
The  averaged  computed  values  are  around  10.5  mm,  11.6  mm,  and  13.2  mm  for  the  three  test 
cases,  respectively,  which  are  comparable  to  the  experimentally  reported  value  of  12  mm. 

Next  the  flame/ wake  dynamics  are  discussed  to  illustrate  the  effect  on  the  Damkohler  number 
similarity.  Only  the  first  two  flames  will  be  compared  here  for  brevity  Similar  conclusions  can 
be  derived  for  the  third  case.  In  the  low  speed/slow  chemistry  flame  (case  1),  the  RZ  tip  is 
found  to  be  oscillating  with  a  10  Hz  frequency  in  the  spanwise  direction.  Figure  15  shows 
eight  instantaneous  snapshots  at  the  plane  Y  =  10  mm  for  case  (1)  over  a  complete  cycle.  The 
temperature  field  is  combined  with  the  contour  of  zero  axial  velocity  to  show  the  extension  of  the 
recirculation  zone.  The  flame  boundaries  are  shown  by  plotting  the  contour  of  the  stoichiometric 
mixture  fraction  (, Z  —  0.021).  As  the  RZ  tip  oscillates,  the  strain  rate  changes  across  the  flame 
surface.  The  RZ  oscillation  combined  with  the  external  shear  layer  (at  the  flame  edge)  result 
in  burned  gasses  shedding,  where  alternative  regions  of  hot  and  fresh  gases  are  separated  at  the 
end  of  the  RZ  and  convect  downstream.  The  unburned  regions  can  grow  in  size  as  they  progress 
downstream  or  they  can  remix  with  hot  gases  and  burn.  The  RZs  in  these  cases  provide  sufficient 
residence  time  to  mix  the  fresh  gases  and  allow  for  burning.  Starting  from  Fig.  15a,  the  RZ 
engulfs  hot  gases  from  its  tip  and  recirculates  hot  products  from  the  left  side  to  burn  the  fresh 
gases  on  the  right  (see  Fig.  15a~15f).  This  in  turn  allows  the  hot  gases  from  the  left  side  to 
convect  downstream  and  mix  with  the  fresh  gas  coming  with  the  shear  layer.  At  the  edge  of  the 
RZ,  quenching  spots  show  up  again  due  to  the  high  strain  rate.  These  quenching  spots  from 
the  left  side  of  the  RZ  shed  downstream  and  mix  with  the  downstream  hot  gases  as  shown  in 
Fig.  15d-15h.  At  the  end  of  the  cycle,  the  RZ  mixes  burned  gases  from  the  right  with  the  fresh 
incoming  gases  from  the  left  (Fig.  15f-15i).  This  process  is  repeated  every  cycle  on  both  sides 
to  act  as  a  natural  self-ignition  mechanism. 

The  same  analysis  is  performed  for  test  case  (2)  in  Fig.  16.  The  flame  shows  wake  and  flame 
structure  similar  to  the  first  flame.  The  RZ  extends  around  50  mm  above  the  bluff-body  base. 
The  oscillation  frequency  is  higher  (about  50  Hz)  as  the  residence  time  is  smaller  in  this  case. 
However,  the  RZ  shows  another  mode  of  instability  where  the  flame  is  detached  from  the  base, 
then  reattaches  at  a  15  Hz  frequency  (Fig.  16g).  This  detachment  is  a  result  of  the  higher  axial 
pressure  gradient  which  overcomes  the  adverse  pressure  gradients  at  every  cycle.  The  flame 
shows  less  coherent  structures  than  test  case  (1).  As  shown  in  Fig.  15  and  Fig.  16,  as  long  as 
the  two  flames  show  similar  stability  dynamics,  they  will  exhibit  Damkohler  number  similarity 
for  the  mean  flame  and  wake  structures. 

Conclusion  -  The  validity  of  Damkohler  number  similarity  for  static  stability  in  augmentor 
flows  is  shown.  The  high-speed/fast  chemistry  flame  shows  similar  recirculation  zone  features 
and  flame  structure  as  the  low-speed/slow  chemistry  flames.  A  characteristic  Damkohler  number 
is  defined  based  on  the  maximum  progress  variable  production  rate.  The  flame  stability  and 
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Figure  15:  Instantaneous  temperature  distributions  for  case  (1)  at  Y  =  10  mm  combined  with 
the  contour  of  zero  axial  velocity  (blue  line)  and  the  contour  of  stoichiometric  mixture  fraction 
(black  line)  over  a  complete  cycle 
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Figure  16:  Instantaneous  temperature  distributions  for  case  (2)  at  Y  =  10  mm  combined  with 
the  contour  of  zero  axial  velocity  (blue  line)  and  the  contour  of  stoichiometric  mixture  fraction 
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main  characteristics  are  found  to  be  controlled  mainly  by  the  wake  dynamics. 

4.2.2  Large-Eddy  Simulations  of  AFRL  Flame  Holder  Experiment 

The  present  section  describes  several  simulations  that  have  been  performed  for  the  planned  ex¬ 
periments  at  AFRL.  The  aim  here  is  to  assess  flow  and  flame  stability  and  to  compare  with 
experimental  data.  Several  simulations  have  been  performed  to  compare  different  flow  and  com¬ 
bustion  effects  and  turbulence  models.  We  have  performed  reactive  and  non-reactive  simulations 
using  a  structured  LES  flow  solver,  where  the  test  rig  geometry  and  the  geometry  of  the  flame 
holder  are  considered  using  the  immersed  boundary  method.  It  is  shown  that,  as  expected  from 
the  flow  regime,  both  the  non-reactive  and  the  reactive  simulations  lead  to  a  stable  flow.  We 
have  further  performed  non-reactive  simulations  on  unstructured  meshes  using  the  CharLES 
flow  solver.  Further,  we  have  performed  unsteady  2D  RANS  simulations  to  assess  the  poten¬ 
tial  of  low  fidelity  simulations.  It  is  found  that  the  2D  URANS  simulations  incorrectly  show 
vortex  shedding  even  for  the  non-reactive  flow.  In  the  following,  we  will  first  describe  the  flow 
conditions  that  have  been  chosen  for  all  simulations  followed  by  the  simulation  results. 

Flow  Conditions  -  Flow-  conditions  have  been  chosen  to  fall  in  a  range  that  is  stable  both 
for  combustion  and  flow.  A  stoichiometric  propane/air  mixture  is  supplied  at  the  inlet.  Inflow 
velocity,  kinematic  viscosity,  and  inlet  temperature  are  specified  as 


^inflow, bulk  *  45.72  TO,/ S 

(25) 

v  -  2,0713  - 10”5  m2/s 

(26) 

^inflow  —  433.0  K  . 

(27) 

These  parameters  lead  to  a  DeZubay  parameter  of  approximately  30. 

Structured  Flow  Solver  Cold  Flow  Results  -  The  structured  mesh  uses  133  x  126  x  108 
cells  in  the  x,  yt  and  z  directions,  and  is  shown  in  Fig.  17.  Figure  18  shows  the  isosurface  of 
the  zero  axial  velocity.  A  two-dimensional  snapshot  of  the  axial  velocity  combined  with  the 
recirculation  zone  iso-surface  is  shown  in  Fig.  19.  The  corrsponding  time  averaged  axial  velcoity 
and  RZ  iso-surface  are  shown  in  Fig.  20,  and  the  mean  profile  of  the  axial  velocity  is  given  in 
Fig.  21. 

Structured  Flow  Solver  Reactive  Flow  Results  -  For  this  case,  the  Inflow  conditions 
are  the  same  as  for  the  non- reactive  case.  A  combined  level  set/flamelet  progress  variable 
model  by  Knudsen  and  Pitsch  [17]  developed  at  Stanford  under  AFOSR  funding  is  used  for 
these  simulations.  The  laminar  flame  speed  for  the  computed  conditions  is  85.933  m/s,  and 
the  laminar  flame  thickness  is  0.398761mm.  Figure  22  shows  an  instantanenous  snapshot  of 
the  temperature  distribution.  Figures  23  and  24  show  instantanenous  snapshots  of  the  carbon- 
monoxide  and  the  progress  variable  distribution,  respectively,  It  is  obvious  from  these  figures 
that  because  of  the  long  residence  time  withing  the  V-gutter  cavity,  the  mixture  is  fully  reacted 
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(a)  Y-Z  plane 


(b)  Z-X  plane 


(c)  3D  cut  view  from  the  outflow  direction 


Figure  17:  AFRL  V-gutter  geometry  and  structured  mesh 


30 


( 


Figure  18:  Instantenous  snapshot  of  the  recirculation  zone  (U  =  Om/s  iso-surface) 


Figure  19:  Instantenous  snapshot  of  the  axial  velocity  component  at  the  Z  =  0.0  plane 
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Figure  20:  Time  averaged  axial  velocity  component  at  the  Z  —  0,0  plane  combined  with  the  RZ 
isosurface 
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Figure  21:  Time  averaged  axial  velocity  profile  at  Z  —  0.0  and  Y  =  0.0 
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Figure  22:  Instantaneous  snapshot  of  temperature  distribution  at  plane  Z  —  0 
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Figure  24:  Progress  variable  at  plane  Z  —  0. 


and  has  therefore  almost  no  CO-  However,  in  the  region  downstream  of  this,  the  finite  reaction 
progress  is  evident  from  the  high  CO  concentrations.  Figure  21  shows  a  comparison  of  the  mean 
axial  velocities  along  the  centerline  from  the  non-reactive  and  the  reactive  simulations*  Both 
velocity  profiles  are  very  similar  up  to  the  end  of  the  V-gutter*  After  that  the  reactive  flow  is 
accelerated  by  the  heat  release. 

CharLES  Low  Mach  Formulation  Cold  Flow  Results  -  The  large  eddy  simulation  (LES) 
code,  CharLES  has  been  used  to  generate  cold  flow  results  for  the  AFRL  flame  holder  geometry. 
An  unstructured,  hexahedral  mesh  of  5.9  million  cells  is  used  for  the  simulation.  Figure  25 
shows  the  computation  domain.  The  underlying  mesh  is  shown  in  Figs.  26,  27,  and  28*  Figures 
29  and  30  show  the  mesh  at  the  Z  —  0  midplane.  Figures  31,  32,  and  33  show  contour  plots  of 
instantaneous  axial  velocity*  time  averaged  axial  velocity,  and  instantaneous  vorticity  magnitude 
in  the  Z  —  0  midplane*  Figure  35  shows  the  vertical  profiles  (y  direction)  of  the  time  averaged 
axial  velocity  at  X— 0.01m  (right  after  the  V-gutfcer),  0.1m  and  0.7m  in  the  Z=0  midplane,  the 
exact  location  being  shown  in  Fig.  34.  It  clearly  shows  the  formation  of  the  wake  in  the  region 
downstream  of  the  V-gutter.  Further  downstream,  the  profile  becomes  more  uniform. 

2D  Unsteady  RANS  Cold  Flow  Results  -  A  leading  commercial  CFD  code  was  used 
to  solve  the  problem  using  an  unsteady  RANS  approach  in  2  dimensions*  An  unstructured, 
hexahedral  mesh  of  around  255,000  ceils,  shown  in  figures  36,  37,  and  38  was  used*  Figures  39 
and  40  show  the  instantaneous  and  the  time  averaged  axial  velocity,  while  Fig.  41  shows  the 
instanteneous  vorticity  magnitude*  These  results  show  coherent  structures  as  the  vortices  are 
shed  behind  the  V-gutter.  This  obviously  is  the  incorrect  solution  in  the  considered  regime. 
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Figure  25:  Computational  domain  for  CharLES  simulation 


Figure  26:  Computational  mesh  for  CharLES  simulation 
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Figure  27:  Computational  mesh  for  CharLES  simulation  (magnified  view) 


Figure  28:  Computational  mesh  near  the  leading  edge  for  CharLES  simulation 
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Figure  29:  Computational  mesh  at  the  Z  =  0  plane  for  CharLES  simulation 


Figure  30:  Computational  mesh  at  the  Z  =  0  plane  for  CharLES  simulation  (magnified  view) 
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Figure  31:  Instantaneous  snapshot  of  the  axial  velocity  component  (m/s)  at  the  Z  —  0.0  plane 
from  CharLES  simulation 
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Figure  32:  Time  averaged  axial  velocity  component  (m/s)  at  the  Z  =  0.0  plane  from  CharLES 
simulation 
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Figure  33:  Instantaneous  snapshot  of  vorticity  magnitude  (1/s)  at  Z  =  0.0  plane  from  CharLES 
simulation 


Figure  34:  Location  of  the  X=0,01m,  QJm  and  0  Jin  lines  in  the  Z  —  0  plane 
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Figure  35:  Profiles  of  time  averaged  axial  velocity  component  (m/s)  at  X=0.01m,  0.1m  and 
0.7m  in  Z  —  0  plane 


Figure  36:  Computational  mesh  for  2D  unsteady  RANS  simulation 
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Figure  37:  Computational  mesh  for  2D  unsteady  RANS  simulation  (magnified  view) 


Figure  38:  Computational  mesh  near  the  leading  edge  for  2D  unsteady  RANS  simulation 
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Figure  40:  Time  averaged  axial  velocity  (m/s)  for  2D  unsteady  RANS  simulation 
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Figure  41:  Instanteneous  vorticity  magnitude  (1/s)  for  2D  unsteady  RANS  simulation 
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